Microscopic description of complex nuclear decay: multimodal fission 
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Our understanding of nuclear fission, a fundamental nuclear decay, is still incomplete due to 
the complexity of the process. In this paper, we describe a study of spontaneous fission using 
the symmetry-unrestricted nuclear density functional theory. Our results show that the observed 
bimodal fission can be explained in terms of pathways in multidimensional collective space corre- 
sponding to different geometries of fission products. We also predict a new phenomenon of trimodal 
spontaneous fission for some rutherfordium, seaborgium, and hassium isotopes. 
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I. INTRODUCTION 

Seventy years ago Joliot Curie and Savitch [ij discov- 
ered that the exposure of uranium to neutrons leads to 
the existence of lanthanum. Following this finding, Hahn 
and Strassmann ^] proved definitively that bombarding 
uranium with neutrons produces alkah earth elements, 
ushering in what has come to be known as the atomic 
age. The term nuclear fission was coined one year later 
by Meitner and Frisch who explained experimental 
results in terms of the division of a heavy nucleus into 
two lighter nuclei. In 1939, Bohr and Wheeler ^] devel- 
oped a theory of fission based on a liquid drop model, 
fnterestingly, their work also contained an estimate of a 
lifetime for fission in the ground state. Soon afterwards, 
Petrzhak and Flerov [11 presented the first experimental 
evidence for such spontaneous fission (SF). 

While early descriptions of fission were based on a 
purely geometrical framework of the nuclear liquid drop 
model [4| (i.e., shape-dependent competition between 
Coulomb and surface energy), it was soon realized [fl 
that the single-particle motion of protons and neutrons 
moving in a self-deforming mean field is crucial for the 
understanding of a range of phenomena such as fission 
half-lives, mass and energy distributions of yields, cross 
sections, and fission isomers 0, IB- macroscopic- 
microscopic method (MMM) proposed by Swiatecki Q 
and Strutinsky et al. [l^, [IH , quantum shell effects are 
added atop the average (or macroscopic) behavior de- 
scribed by the liquid drop, and this approach turned out 
to be very successful in explaining many features of SF 

Quantum mechanically, fission represents a time- 
dependent solution of the many-body Schrodinger equa- 
tion where all particles move collectively. To fully solve 
such a time-dependent problem for more than 200 parti- 
cles is neither possible nor sensible because the essence of 
the process is in its coherence. Consequently, most of the 
essential physics should be contained in underlying mean 



fields. This determines the choice of a microscopic tool 
to be used: the nuclear density functional theory (DFT) 
p^ . The advantage of DFT is that, while treating the 
nucleus as a many-body system of fermions, it provides 
an avenue for identifying the essential collective degrees 
of freedom. 

Because the commonly used nuclear density function- 
als are usually adjusted to nuclear ground-state proper- 
ties and infinite nuclear matter, and most applications 
are symmetry-restricted to speed up computations, self- 
consistent theories typically are not as quantitative as 
MMM when it comes to SF, except for some cases pjSj . 
It is only recently that an effort has been made to sys- 
tematically optimize the effective forces by considering 
experimental data relevant to large deformations 

In this work, we provide a microscopic description of 
multi-mode fission based on the symmetry-unrestricted 
nuclear DFT. Since many observed fission characteristics 
can be traced back to topologies of fission pathways in 
multidimensional collective space 12], allowing for arbi- 
trary shapes on the way to fission is the key. To this 
end, we search for the optimum collective trajectory in a 
multidimensional space. The barrier penetration proba- 
bility, or a fission half-life, is computed by integrating the 
action along this optimum path. In practice, this is done 
by constraining the nuclear collective coordinates associ- 
ated with shape deformations to have prescribed values 
of the lowest multipole moments, Qx^i by which we ex- 
plore the main degrees of freedom related to elongation 
(A/i — 20), reflection-asymmetry (A/x = 30), and necking 
(A/i — 40). The effects due to triaxiality arc known to 
be important around the top of the first fission barrier 
[H [il [H m [U m. Indeed, the first saddle point is 
lowered by several MeV by triaxial degrees of freedom. 
In our work, we take into account the impact of non-axial 
degrees of freedom by considering the triaxial quadrupole 
moment (A/J, = 22). 

The paper is organized as follows. Section HIl describes 
the details of our model. The results of calculations are 



2 



discussed in Sec. IIIIl Finally, Sec. IIVI summarizes the 
main results of our study and offer perspectives for future 
work. 



II. THE MODEL 

The calculations were carried out using a symmetry- 
unrestricted DFT program based on the Hartree Fock- 
Bogohubov solver HFODD [H ^ capable of treating 
simultaneously all the possible collective degrees of free- 
dom that might appear on the way to fission. Based 
on this DFT framework, we calculated the total energy 
along the fission pathways, corresponding collective iner- 
tia (collective masses) and zero point energy (ZPE) cor- 
rections to account for quantum fluctuations. 

In the particle-hole channel, we use the SkM* energy 
density functional [2^ that has been optimized at large 
deformations; hence it is often used for fission barrier 
predictions. In the pairing channel, we adopted a se- 
niority pairing force with the strength parameters fit- 
ted to reproduce the experimental gaps in ^^^Fm [2^ . 
Because the nuclei considered are all well bound, pair- 
ing could be treated within the BCS approximation. 
The single-particle basis consisted of the lowest 1,140 
stretched states originating from the lowest 31 major os- 
cillator shells. 

In the analysis of fission pathways, we explored mul- 
tidimensional collective space. To separate fission path- 
ways, we computed energy surfaces in the deformation 
spaces {Q20, Q30} and {Q20, Qio}- The calculations were 
not limited to axial shapes; triaxial deformations appear 
if energetically favorable (e.g., within the inner barrier). 
At each deformation point, defined by the set of con- 
strained multipole moments, fully self-consistent DFT 
equations have been solved, whereupon the total energy 
of the system is always minimized with respect to all 
remaining (i.e., unconstrained) shape parameters. The 
optimum ID paths have been localized in the form of 
multipole moments, Q30, Q22, and Qaq, becoming func- 
tions of the driving moment, (520 • 

The vibrational and rotational ZPE corrections and 
the cranking quadrupole mass parameter were calculated 
as described in Ref. [l^]. The spontaneous fission half- 
lives were estimated from the WKB expression for the 
double-humped potential barrier [2^ [2^ assuming a ID 
tunneling path along Q2q- 

III. RESULTS 

To demonstrate the validity and generality of our 
method, we chose a case where several fission pathways 
were known to coexist and all intrinsic symmetries of the 
nuclear mean field were broken. In this respect, a phe- 
nomenon known as bimodal fission, observed in several 
fermium and transfermium nuclei |30l . Isil [3^ . Issj . is a 
perfect testing ground. It manifests itself, for example. 



in a sharp transition from an asymmetric mass division 
in ^^^Fm and ^^^No to a symmetric mass split in ^^*Fm 
and ^^^No. Furthermore, the total kinetic energy dis- 
tributions of the fission fragments appear to have two 
peaks centered around 200 and 233 MeV. It has been 
suggested [H, HI IH, [H, [Ullil that the higher energy 
fission mode corresponds to a scission configuration as- 
sociated with two touching, nearly spherical, fragments 
with the maximal Coulomb repulsion, whereas the lower- 
energy mode can be associated with more elongated frag- 
ments. Before this work, bimodal fission was studied 
within the MMM 111, H HI, [H, [13 and nuclear DFT 
[ll[23,[3i,[43,lii!,[43. All those studies were symmetry- 
restricted (i.e., they did not consider simultaneous inclu- 
sion of elongation, triaxiality, and reflection-asymmetry). 



A. Fission pathways of the fermium isotopes 




FIG. 1: (Color online) Two-dimensional total energy surfaces 
for ^^*Fm in the plane of collective coordinates: Q20-Q30 (a), 
and Q20-Q40 (b). The fission pathways are marked: symmet- 
ric compact fragments (sCF), symmetric elongated fragment 
(sEF), and asymmetric elongated fragments (aEF) pathways. 
The difference between contour lines is 5 MeV in (a) and 
10 MeV in (b). The asymmetric trajectory aEF bifurcates 
away near Q20 ~ 150 b from Q3o=0 while the bifurcation be- 
tween sEF and sCF occurs near Q20 ~ 225 b. The inset shows 
the multipole moments Qxo along sCF and sEF. 
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To identify saddle points and fission pathways in a mul- 
tidimensional energy surface is not a simple task. As 
pointed out in earlier studies [H, [13, Hj, [31], saddle 
points obtained in calculations constrained by only one 
collective variable are sometimes incorrect; hence, special 
numerical techniques are required to find them. To this 
end, we computed two-dimensional (2D) energy surfaces 
in {Q20, Q30} and {Q2Q, Qao} planes. Based on the ini- 
tial 2D uniform grids, Q20=0(10)400b, Q30 =0(5)50b3/2, 
and (540 = 0(10)200 b^, we calculated the constrained 
HF-I-BCS energy. In our HFODD calculations we em- 
ployed the standard method of q uadratic constraint (the 
quadratic penalty method) [U [4^. Within this ap- 
proach, the moments calculated from the converged den- 
sities differ slightly from the requested values defining the 
constraints [46] . This implies that the final mesh used for 
interpolation is non-uniform. Using the self-consistent 
values of multipole moments, the non-uniform interpola- 
tion was carried out to produce the final result. Figure [1] 
shows 2D energy surfaces for ^^^Fm. 

The ID fission pathways shown in Fig. [T] were obtained 
by finding the local energy minimum at a given value of 
Q20 and following this solution when gradually increas- 
ing the driving moment Q2o- This method thus relies on 
local properties of the (520-constrained energy surface. 
In practice, one obtains different energy surfaces locally 
valid around each fission pathway. In order to facilitate 
the presentation, however, we carried out interpolation 
between different energy sheets. For that reason, the con- 
tour lines in Fig. [T] slightly depend on the algorithm used 
to interpolate between 2D constrained results; hence, the 
2D surfaces should be considered as a qualitative guid- 
ance on the topological structure of the energy surface. 

Beyond the first barrier, at Q20 ~ 150 b, a refiection- 
asymmetric path corresponding to asymmetric elongated 
fragments (aEFs) branches away from the symmetric val- 
ley, see Fig. [T^. At Q20 — 225 b, a refiection-symmetric 
path splits into two branches: one corresponding to 
nearly spherical symmetric compact fragments (sCFs) 
and one associated with symmetric elongated fragments 
(sEFs). This bifurcation is clearly seen in Fig. [TJj. Such 
three fission pathways were predicted in early work based 
on MMM [34,131] and also found recently within a DFT 
framework [40j, 1411] , except that the axial-symmetry was 
enforced in all these studies. 

A pattern of similarly competing fission valleys was 
found for all investigated isotopes. As an illustrative 
example, the results of our ID fission pathway calcu- 
lations for 252,258,264pj^ displayed in Fig. [1 (The 
energy curves corresponding to individual pathways in 
256,260pj-|^ can be found in our earlier work 26].) The 
transition from an asymmetric fission path in ^^^Fm to a 
compact-symmetric path in ^^^Fm is due to shell effects 
in the emerging fission fragments approaching the dou- 
bly magic nucleus "^Sn [ll [H, [13, (3, [!§] . The triaxial 
deformations are important around the first (inner) fis- 
sion barrier, and they reduce the fission barrier height by 
several MeV. 
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FIG. 2: (Color online) Calculated ID fission pathways for 
252, 258, 264 Pjj^ as functlons of the driving quadrupole moment, 
Q20- Open circles, light diamonds, and dark diamonds denote 
the symmetric compact fragment (sCF) , symmetric elongated 
fragment (sEF), and asymmetric elongated fragment (aEF) 
valleys, respectively. The effect of triaxiality is important in 
the region of the first barrier; the corresponding energy gain 
is marked by gray shading. 



While the asymmetric pathway aEF is favored in the 
lighter Fm isotopes, e.g., ^^^Fm, both symmetric paths 
are open for ^^'^Fm, due to the disappearance of the outer 
fission barrier in sCF and sEF. It is to be noted that the 
symmetric pathways sCF and sEF in ^^^Fm are predicted 
to bifurcate away well outside the first barrier (see also 
recent work [FO] based on MMM). 

In the case of 260-264-pj-|^^ g^^j ^j-^^^^ there is no outer 
potential barrier along the sCF trajectory, and the sEF 
and aEF paths lie significantly higher in the outer region. 
It should be emphasized that the pathways correspond to 
different regions of the collective space and this is appar- 
ent when studying them in more than one dimension. 
Indeed, aEF is well separated from sCF and sEF in Q30 
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(the apparent crossing between aEF and sEF in Fig. [T]3 
is an artifact of the 2D projection) while the symmetric 
trajectories sCF and sEF strongly differ in the values of 
higher multipole moments A=4, 6, and 8 (see the inset 
in Fig. [J). 



B. Spontaneous fission half-lives 
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FIG. 3: (Color online) Potential energy curves of ^^^Fm, 
^^®Fm, ^^®Fm, and ^''"Fm along optimum ID paths contain- 
ing the quantum zero-point energy correction, drawn in the 
common scale relative to values calculated at the ground-state 
minima. 



Having determined the lowest fission valleys, i.e., the 
lowest energy pathways along Q20j we computed the cor- 
responding ZPE corrections. The resulting collective po- 
tentials are plotted in Fig. [3] for ^^^Fm, ^sepm, ^^^Fm, 
and ^^°Fm. It is seen that when going from ^^^Fm to 
260pj^^ the first barrier gets reduced and the outer bar- 
rier disappears altogether. 

To assess the SF half-lives theoretically, we calculated 
the collective inertia parameter along Q20 and performed 
WKB barrier penetration calculations for even-even fer- 
mium isotopes with 242<A<260. We assumed two val- 
ues of the ground state energy counted from the ground 
state potential energy minimum: £^0=0.3 MeV and the 
commonly used value [sH of 0.5 MeV. The resulting SF 
half-lives are shown in Fig. 3] In spite of a fairly simple 
ID penetration picture, it is satisfying to see a quantita- 
tive agreement between experiment |52l . [ssj and theory 
(for Gogny-DFT results, see [i^l)- The existence of a 
small outer barrier in ^^^Fm is significant as it increases 
the fission half-life in this nucleus by more than four or- 
ders of magnitude compared to that of ^^^Fm, thus ex- 
plaining the rapid change in experimental SF half-lives 
between these nuclei. 
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FIG. 4: (Color online) Fission half lives of even-even fermium 
isotopes with 242<A<260, calculated in this study compared 
with experimental data j52|, for the two values of the zero- 
point energy Eo^O.SMeV and 0.5 MeV. 



C. Multimodal fission 

To map out the competition between different fission 
pathways in the heaviest elements, we carried out sys- 
tematic calculations for even-even nuclei with 98<Z<108 
and 154<Af<160. A transition from the usual asymmet- 
ric fission channel seen in the actinides to compact sym- 
metric fission is seen when moving towards ^^"^Fm. In 
the intermediate region of bimodal fission, two symmet- 
ric channels coexist. Around ^^''Sg (Z=106, A^=154), our 
calculations predict trimodal fission, i.e., competition be- 
tween the asymmetric fission valley and two symmetric 
ones. (The term "multimodal fission" has been previ- 
ously used by M.G. Itkis et al. [54] in the context of 
fusion-fission and quasi-fission of hot superheavy nuclei 
produced in heavy ion reactions.) 

The representative fission pathways for ^^"'Cf (asym- 
metric fission), ^^"^Fm (symmetric compact fission), 
258]\fQ (bimodal fission) and ^^^Hs (trimodal fission) are 
displayed in Fig. O The inclusion of triaxiality signif- 
icantly reduces the inner barrier in ^^''Cf, ^^°Fm, and 
258]\fQ "vvhile the effect in ^^^Hs is much weaker. Other 
nuclei from this region exhibit similar pattern of fission 
pathway competition. A summary of our findings is given 
in Fig. [51 



IV. SUMMARY 

In summary, the symmetry-unrestricted nuclear DFT 
framework has been applied to the problem of SF. As 
an example, we studied the challenging case of static 
SF pathways in ^^^Fm, ^^^Fm, and ^^°Fm and in a 
number of neighboring nuclei. We found competition 
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FIG. 5: (Color online) Similar as in Fig. [2] except for fission 
pathways in ^^*Cf (asymmetric fission), ^^°Fm (symmetric 
compact fission), ^^*No (bimodal fission) and '^^'^Hs (trimodal 
fission) . 



between symmetric-compact, symmetric-elongated, and 
asymmetric-elongated fission valleys that is consistent 
with the observed distribution of fission yields. The sad- 
dle points obtained in constrained ID calculations were 
confirmed through an analysis of 2D energy surfaces. 
From the calculated collective potential and collective 
mass, we estimated SF half-lives, and good agreement 
with experimental data was found. Finally, we predicted 
trimodal fission for several rutherfordium, seaborgium, 
and hassium isotopes. 

It is worth noting that calculations of self-consistent 
energy 2D surfaces are computer intensive. Because a 
single HFODD run with all self-consistent symmetries 
broken takes about 60 minutes of CPU time, it takes 
about 3 CPU-years to carry out the full fission pathway 
analysis for 24 nuclei; hence, massively parallel computer 
platforms had to be used. 

In the near term, we intend to improve the theory of 
SF half-lives by considering multidimensional inertia ten- 
sors and by performing the direct minimization of the 
collective action in a multidimensional collective space 
[sst - In the long term, the theory will be extended to 
account for nonadiabatic effects (e.g., along the lines of 
Refs. [sgUst}). In addition, quality microscopic input for 
fission calculations is needed. Of particular importance 
is the development of the nuclear energy density func- 
tional better reproducing both bulk nuclear properties 
and spectroscopic data. 
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